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r"| ■ Abstract. Currently used model of spherical accretion onto a magnetized rotating neutron star 

fS I encounters major difficulties in explaining the entry rate of accreting material into the stellar field 

and spin evolution of long-period X-ray pulsars. These difficulties can be, however, avoided if 
the magnetic field of the material captured by the neutron star is incorporated into the model. 
The magnetic field of the flow itself under certain conditions controls the accretion process and 
significantly affects the parameters of the accreting material. The mode by which the accretion flow 
enters the stellar magnetosphere in that case can be associated with Bohm (or turbulent) diffusion 
and the torque applied to the neutron star appears to be substantially higher than that evaluated in 
the non-magnetized accretion scenario. 
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)Q : 1. INTRODUCTION 

o 

High Mass X-ray Binaries (HMXBs) are interacting binary systems made of a massive 
O/B-type star and a magnetized, rotating neutron star. The massive star underfills its 
Roche lobe and looses material in a form of stellar wind. As the neutron star moves 
through the wind of density with a relative velocity V K \ = V ns + V W , it interacts in a 
unit time with the mass 9jt c = 7lR^pJV m \. Here R ns and M ns are the radius and the mass 
of the neutron star, V ns is the velocity of its orbital motion, and V w the wind velocity. The 
Bondi radius, Rq = 2GM ns /V^, represents the maximum distance at which the neutron 
star is able to capture material from its environment. 

The captured material initially follows ballistic trajectories and approaches the neu- 
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tron star to the distance Ra = ^/i 2 /9^c(2GM ns ) 1//2 j at which the magnetic pressure 

due to the dipole field of the star reaches the ram pressure of the accreting material. Ra 
represents the Alfven radius and \x is the dipole magnetic moment of the neutron star. 
The accretion process inside the Alfven radius is fully controlled by the stellar field. The 
accreting material penetrates into the field at the magnetospheric boundary and reaches 
the stellar surface flowing along the magnetic field lines to the magnetic pole regions. 



TABLE 1. Persistent Long-period X-ray Pulsars* 



Name Sp. type 1 " P orb , d P s , s P, s/s * CRSF logL x i d 



GX 301-2 


Bl la 


41.5 


683 


[±] 


E- 


-7.3 30-38keV 


37- 


37.5 


1.8-3 kpc 


X Persei 


BO Ve 


250 


837 


[±] 


E- 


-5.5 29keV 


34.7 


-35.5 


950 pc 


4U 2206+54 


09e 


19.25 


5554 


[ + ] 


E- 


-6.3 


35- 


35.6 


2.6 kpc 


2S 01 14+65 


Bl la 


11.6 


10008 


[-] 


E- 


-6.1 


35 


-36 


3 kpc 



* For references see [18] 

t Spectral type of the massive star, 

$ Spin behavior ( [ — ] for spin-up and [ + ] for spin-down) and maximum absolute value of P, 
H L x is in ergs -1 



The accretion power is then converted into X-rays. If the star is obliquely rotating it ap- 
pears as an accretion-powered pulsar. The period of pulsations in this case is equal to the 
spin period of the neutron star. The mass capture rate by the neutron star in a stationary 
accretion picture can be evaluated through observations of the source X-ray luminosity, 
L x , 2LsWl a = L x R RS /GM ns . 

The accretion flow beyond the magneto spheric boundary can be treated in spherically 
symmetrical or disk approximation. The specific angular momentum of the captured ma- 
terial is fc a o ~ Q. 0I \,Rq, where ^ or b = 2n/P 0I h is the orbital angular velocity and ZVb the 
orbital period of the binary system. The transfer rate of angular momentum associated 
with accretion is J = ^Q.^RqTI^ where ^ is the factor by which angular momentum 
accretion rate is reduced due to inhomogeneities (velocity and density gradients) and 
magnetic viscosity in the accretion flow. The angular velocity of the material increases 
during accretion process from the initial value C0q = Q or b to 0) en (r) = ^0)o (Rq/t) 2 . As 
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it reaches the Keplerian angular velocity, 00^ = (GM ns /r 3 ) , the accretion disk forms. 
The distance at which formation of the disk in a wind-fed binary system can be expected 
is referred to as the circularization radius, 7? c j rc = J 2 / GM ns 9jt 2 . It is defined by equating 
Wen(^circ) = Wk(^circ)- The Keplerian disk around a magnetized neutron star could be 
formed only if R c i IC > Ra- Otherwise, the stellar magnetic field would prevent the ma- 
terial from reaching the circularization radius. This inequality can be expressed as (see 
[18], and references therein) 

y rel < V cr c 200 kms- 1 $J/ 4 ^J/ 14 m 11 / 28 P^ 4 w\ /2 \ 

where /I30 = ju/ 10 30 Gem 3 and m is the mass of the neutron star in units of 1 .4Mq. The 
system orbital period, P40, is expressed in 40 days, and Win = 9?lc/10 17 gs _1 . Finally, 
<d0.2 = <^/0-2 is normalized to its maximum average value in the case of non -magnetized 
accretion flow (see e.g. [33], and references therein). 

The wind velocity of O/B-type stars usually ranges between 400 to 1000 kms -1 . A 
smaller wind velocity can be expected only in Be-type stars, which are surrounded by a 
viscous outflowing disk. But even in that case the orbital velocity of neutron stars would 



exceed V CI if the orbital period of the binary system satisfies the condition P or b < 50 days. 
It, therefore, appears that a formation of Keplerian accretion disk in a large number of 
HMXBs is not expected and the accretion process in those systems can be considered in 
a quasi-spherical approximation. 

About 60 X-ray pulsars in our Galaxy are associated with HMXBs [24]. More than 10 
of them represent the subclass of persistent long-period X-ray pulsars. Parameters of the 
best studied members of this subclass are summarized in Table 1 . They are the wind-fed 
accretors. Observations give no evidence for a presence of persistent accretion disk in 
these systems. Instead, there are some indications that the neutron star in 2S 01 14+65 
is embedded in and accreting from a dense low angular momentum gas cloud [27]. 
GX 301-2 is the bright pulsar with the X-ray luminosity (1 -3) x 10 37 erg s _1 . The rest 
of the sources are fainter: L x < 10 36 ergs _1 . Nevertheless, the X-ray spectra of all these 
pulsars look similar, which strongly suggests that the accretion process operating in 
these binary systems is governed by the same mechanism. Observations of the cyclotron 
resonance scattering feature (CRSF) in X-ray spectra of GX 301-2 and X Persei suggest 
that the surface field of the neutron stars lies in the range of (2 — 4) x 10 12 G. All of these 
pulsars display a spin evolution. 2S 01 14+65 shows spin-up and 4U 2206+54 spin-down 
behavior. Both the spin-up and spin-down phases have been observed in GX 301-2 and 
X Persei. The recent spin history of GX 301-2 is shown in Fig 1. The long-term spin- 
down trends at the rate V s d ~ — 10 _13 Hzs _1 are superposed on slow spin-up trends and 
rapid spin-up episodes at the rate v su ~ 5 x 10 _12 Hzs _1 [25], [20]. 

The accretion picture in the pulsars described above is the subject of our paper. 
We show that the currently used non-magnetized accretion flow scenarios encounter 
major difficulties in explaining the spin evolution and persistent character of the pulsars 
(next Section). Then we discuss the strength of the magnetic field in the accretion 
flow, elaborate the Magnetically Controlled Accretion (MCA) scenario and consider 
its application to the interpretation of GX 301-2. A brief summary of our results is 
concluding the paper. 



The magnetic field of the accretion flow is usually neglected in currently used spherical 
accretion scenarios. It is assumed that the neutron star captures material at the rate 
Wl a = Tl c , which is then falling towards the star in spherically symmetric fashion 
with the free-fall velocity Vff(r) = (2GM m /r) 1 ^ 2 . The gravitational energy of the star 
is converted predominantly into kinetic energy of the accreting material and its ram 
pressure increases as 



where £ r o = P«y%>\ is the ram pressure at the Bondi radius. 

The flow is decelerated by the stellar magnetic field when it approaches the star to a 
distance of ~ Ra- This leads to the formation of a magnetosphere of radius R m and to the 
heating of the flow up to a temperature T s = (3 /\6)Tff(R m ), where Tff(r) = GM ni m v /k^r 
is the proton free-fall temperature, m p is the proton mass and k% is the Boltzmann con- 



2. NON-MAGNETIZED ACCRETION FLOW SCENARIO 




(1) 




FIGURE 1. Pulse frequency and 12-50 keV pulsed flux of GX 301-2 observed with the Fermi Gamma- 
ray Burst Monitor detectors. GMB Pulsar Project http://gammaray.nsstc.nasa.gov/gbm/science/pulsars 



stant [22]. The heating occurs on the dynamical (free-fall) time, tff(r) = (r 3 /2GM ns ) 
in the shock located at the magneto spheric boundary. The magnetosphere surrounded 
by the hot spherically symmetric accretion flow is closed (i.e. does not contain open 
field lines). It is rotationally symmetric about the magnetic axis and reflection symmet- 
ric about the equatorial plane. The boundary of the magnetosphere is convex towards the 
accreting material and contains two cusp points situated on the stellar magnetic axis. The 
size of the magnetosphere in the plane of the magnetic equator is R m ~ (1.7 — 2) x 7? a 
and the distance to the cusp points is close to the Alfven radius of the neutron star (see 
[1], [11], [29]). 

Formation of the magnetosphere, in the first approximation, prevents the accret- 
ing material from reaching the stellar surface. To enter the magnetosphere the ma- 
terial should move across the magnetic field lines. However, the rate of diffusion, 
M m ~ 2^7? m 5 m poVff (see e.g. [12], [19]), in the considered case is limited to M- m ~ 
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(<5m/7? m )9JI c , where 5m = (D e ff z^) lS the thickness of the diffusion layer, D eff is the 
effective diffusion coefficient, Td is the diffusion time and po = Wl c / (AnR^Vff) is the 
density of the free-falling material at the magneto spheric boundary. In particular, the 
rate of Bohm diffusion, D e ff ~ Db = O^ck^T^/eB, on the dynamical timescale, Td = tff, 
is limited to 

M B < 6 x K^gs- 1 x c4 /2 m 1 / 7 ^ 14 . (2) 

which is a few orders of magnitude smaller than the mass transfer rate expected in the 
case of stationary accretion picture [17]. 



A higher entry rate could be expected if the boundary were interchange unstable. The 
fastest mode is the Rayleigh-Taylor instability (see [1], [2], [11], [3]). However, this 
instability in the considered case can occur only if the temperature of the material at the 
boundary is limited to T(R m ) < 0.1 — 0.3Tff(R m ) [1], which implies that cooling of the 
accreting material dominates heating. Otherwise, the effective acceleration applied to 

GM ns V% 

the material at the boundary, g e ff = 2 cos^ — — l —, will be directed outwards 

from the star. Here 6 is the magnetic latitude, % is the angle between the radius-vector 
and the normal to the magneto spheric boundary, and R curw is the curvature radius of the 
field lines. Studies [11], [3] have shown this criterium to be satisfied only in bright X-ray 
pulsars with L x > L cr , where 

L cr ~ 3 x \& 6 ^m l ' 2 R~ 1/8 ergs" 1 , (3) 

and 7?6 = R ns /10 6 cm. Otherwise, cooling time of the material at the magneto spheric 
boundary (which is governed by the inverse Compton scattering of hot electrons on X- 
ray photons emitted from the surface of the neutron star) is larger than the dynamical 
time on which the heating of the electrons occurs. The faint X-ray pulsars within this 
scenario can operate only as transient burst-like sources [22]. 

The Kelvin-Helmholtz instability can be effective in fast spinning pulsars with P s ~ 

3 1 mT 5 / 1 9# 15 3 ^ 7 s, in which the radius of magnetosphere is close to a so called co- 
rotation radius of the neutron star, R COI = (GM ns / (Os) 1 / 3 , at which the linear velocity at 
the magneto spheric boundary is equal to the Keplerian velocity [8]. Finally, the turbulent 
diffusion can also be effective only in fast spinning pulsars since the velocity of the 
turbulent motions is limited to the relative velocity between the accreting material and 
magnetosphere, which does not exceed ~ C0 s R m [9]. 

Thus, the currently used spherical accretion scenario onto a magnetized neutron star 
predicts that faint long-period X-ray pulsars should be transient burst-like sources in 
which the accretion process onto the stellar surface occurs on the cooling time of the 
material at the magneto spheric boundary. This prediction contradicts observations. 

The spin evolution of an accretion-powered pulsar is described by the equation 

I(Q S = K SU -K sd , (4) 

where / is the moment of inertia of the neutron star. The spin-up torque applied to the 
star from the quasi-spherical accretion flow is K su = ^k a oTl a = ^Q. 0I \,RQWl a . The spin- 
down torque can be evaluated by solving a model task in which the magnetosphere is 
approximated by a sphere of radius R m and its interaction with surrounding material is 
treated in terms of rotation in viscous medium. The spin-down torque in this case is 
K s d ~ AnR^VtpoVtp . Assuming the turbulent nature of the viscosity, v t = £tR m V(j,, where 
V(j) < C0 s R m is the relative azimuthal velocity between the magnetosphere and material at 
the boundary, one finds K sd = e t /i 2 /R\ 0T [25]. Here £ t is the efficiency parameter. Using 
parameters of GX 301-2 one finds, that the spin-down torque applied to the neutron 
star is a factor of 1000 (B/Bqrsf) smaller than the spin-down torque evaluated from 
observations. Here 5crsf = 4 x 10 12 G is the surface field on the neutron star in this 



system measured through observations of CRSF. This inconsistency indicates that either 
the surface field of the star is a factor of 30 stronger or the accretion picture in this source 
differs from that expected in the currently used spherical accretion scenario. 

3. MAGNETIZED VS. NON-MAGNETIZED ACCRETION FLOW 

Spectropolarimetric observations [14, 31, 26] have shown that a relatively strong mag- 
netization of O/B-type stars is not unusual. The strength of the large-scale field at the 
surface of these objects has been measured in the range B* ~ 100 — 1000 G, and in some 
cases beyond 10 kG. The field in the environments of these objects can be evaluated by 
taking into account that the dipole approximation (B °< r 3 ) to the magnetic field of mas- 
sive hot stars remains valid up to a distance R^, at which the dynamical pressure of their 
stellar wind reaches the magnetic tension. The field in the wind propagating beyond this 
point decreases as B r~ 2 [35]. Therefore, the magnetic energy density in the wind at a 
distance a > 7? k is £ m (a) = 111/ (2nR^a 4 ) , which implies 

£ m (a) ~ 0.03 erg cm" 3 ( —^-) * ( - ^* , ) ' ( rr^—) ' , (5) 
v; 5 V10 13 cm/ V 3xl ° 38 Gcm 3 y \100R & J 

where /i* = (\/2)B*R\ is the dipole magnetic moment of the massive star. 

The thermal energy density of the wind at the Bondi radius, £ t h = p^V 2 , can be 
evaluated by taking into account that p^ = 9Jl c / kR^V k \ as 

£th = 2 ~ 0.02 ergcm- 3 m- 2 Wl l7 ( c .J rel , ) * ( . ) ' . (6) 

7t(GM ns ) 2 s 5 V 500kms / V 106cms_ J 

Here V s is the sound speed in the wind at Rq. Thus, if the neutron star is situated in a 
relatively slow wind the ratio £ t h(^G)/£m(^G) = P can be about or even smaller than 
unity. The magnetic field of the accretion flow cannot in this case be neglected. 

4. MAGNETICALLY CONTROLLED ACCRETION SCENARIO 

We consider spherical accretion onto a magnetized neutron star under the assumption 
that the parameter /3 in the material captured by the star at the Bondi radius is close 
to unity. The ram pressure of the captured material at the Bondi radius significantly 
exceeds its thermal and magnetic pressure. The Alfven velocity, Va = Bf/(Anp) 1 / 2 , 
in the accreting material under these conditions is substantially smaller than the free- 
fall velocity, and, hence, the field annihilation time exceeds the dynamical time. 
This indicates that accretion process can initially be considered in the magnetic flux 
conservation approximation. Here Bf is the field strength in the accretion flow. 

The magnetic field in the free-falling material is dominated by the radial compo- 
nent [36], which under the condition of the magnetic flux conservation increases as 
B r (R) ~ Bf(Rc) (r/7?o) [4]. The magnetic pressure in the accreting material, there- 



fore, increases as it is approaching the star as 



£m(r) = £ m (* G ) (J^-J , (7) 

while the ram pressure of the free-falling spherical flow is £ r oc r ~ 5 / 2 ( see Eq. 1), and 
thus, £ m /£ r oe r~ 3 / 2 . This indicates that the gravitational energy of the neutron star in 
the process of spherical accretion is converted predominantly into the magnetic energy 
of the accreting material. 

The distance 7? s h, at which the magnetic pressure in the accretion flow reaches the ram 
pressure, can be evaluated by equating £ m (^sh) = £r(^sh)- This yields [34], 

WM 4/ % = «-2/3 2GM m V* 13 



rel 

For i? s h (hereafter Shvartsman radius) to exceed Alfven radius of the neutron star, the 
relative velocity should satisfy condition V re \ < V mcSL , where 

y mca = r l/5 ^l /3 \2GM ns yy^ /5 ^ 5 . ( 9) 

For typical parameters of HMXBs one finds 

V mca ~ 680 J8-V5 W lV3 5jU -6/35 ^3/35 y 2/5 

which under the conditions of interest substantially exceeds the value of V CI . Thus, 
one can distinguish a subclass of HMXBs in which the accretion occurs in spherically 
symmetrical fashion and the accreting material is strongly affected by the magnetic field 
of the flow itself. This subclass is defined by condition V CI < V ie \ < V mca . 

Rapid amplification of the magnetic field in the spherical flow as well as deceleration 
of the flow by its own magnetic field at the Shvartsman radius have been confirmed 
by the results of numerical studies of magnetized spherical accretion onto a black hole 
[15, 16]. These calculations have shown that the magnetized flow is shock-heated at the 
Shvartsman radius up to the adiabatic temperature. If the cooling of the material at R s h is 
inefficient the flow switches into the convective-dominated stage in which some material 
is leaving the system in a form of jets and the mass accretion rate in this region proves 
to be smaller than its initial value. 

Another important effect, the Compton cooling of the accreted material, can modify 
this accretion picture. Due to a large value of V tt \, the Shvartsman radius in these systems 
is located much closer to the neutron star which is a powerful source of X-ray emission. 
The Compton cooling of the accreting material in this case can prevent the flow from 
switching into the convective-dominated stage and being ejected out from the binary 
system. 

The plasma cooling at can be effective if the Compton cooling time [11], 



is smaller than the characteristic time of plasma radial motion inside the Shvartsman 
radius. Here m e is the electron mass and c T is the Thomson cross-section. If the magnetic 
flux in the accreting material is conserved, the accretion ends at the Shvartsman radius. 
Further accretion in this case is impossible. Otherwise, the magnetic energy in the 
flow would exceed its gravitational energy, which contradicts the energy conservation 
law (for discussion see [34]). Therefore, the accretion flow can approach the star to a 
closer distance only if a dissipation of the magnetic field in the flow starts. If the field 
dissipation is governed by magnetic reconnection the characteristic time of the accretion 
process inside is limited to t > t iec , where 



-i, (Ya 



? rec = ZTTT = *ff ~ • (U) 



Theoretical studies [32] and observations [30] suggest the value of the efficiency pa- 
rameter r\ m to be in the range of 0.01-0.15. It should be also mentioned that the Alfven 
velocity in the accretion flow reaches the free-fall velocity at the Shvartsman radius. 
Combining Eqs. (10) and (11), one can express inequality t c < t KC as 



L x > ri m L CI ^ j . (12) 

This indicates that the cooling of magnetized flow can be effective even in faint X-ray 
pulsars in which the Shvartsman radius does not significantly exceed the Alfven radius 
of the neutron star. 

The accretion picture of a cold magnetized gas has been discussed in [5,6]. It has been 
shown that the material in this case tends to flow along the field lines and in the region 
r < ^sh accumulates in a dense non-Keplerian slab (see Fig. 1 in [6]). The material in the 
slab is confined by the magnetic field of the flow itself and its radial motion continues 
as the field is annihilating. The accretion process in the slab, therefore, occurs on the 
timescale of t iec . This indicates that the thermal pressure in the slab is equal to magnetic 
pressure and hence, the gas density at the inner radius of the slab, which is equal to the 
magneto spheric radius, can be evaluated as 

/i 2 m p 

Psl= w^- (13) 

Here Tq is the gas temperature at the inner radius of the slab. The thickness of the 
slab depends on the magnetic field configuration and in the first approximation can be 
evaluated from the continuity equation as 

hz(R m ) = ^. (14) 

The rate of plasma diffusion from the slab into the stellar field can be evaluated as 

9Kin(tfm) = 47Ttf m S m p sl V ff = AtzR^D^ p sl y f f /2 , (15) 
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where 5m = (D e ff Td) 1 and the diffusion time, Td = tff(R m ) = R m /Vff(R m ), is equal to 
the dynamical time on which the material inside the magnetosphere leaves the diffusion 
layer by flowing with the free-fall velocity along the stellar field lines towards the stellar 
surface. As long as Wli n is smaller than the mass capture rate and, correspondingly, 
the rate of mass transfer in the slab in the radial direction, Wl a , material accumulates 
at the inner radius of the slab and the gas pressure in this region increases. The inner 
radius of the slab is defined by equating the gas pressure to the magnetic pressure by 
the dipole field of the neutron star. Therefore, as the material at the slab accumulates, 
its inner radius decreases and the material is getting closer to the star. A stationary 
accretion picture is expected as the magneto spheric radius reaches a value at which 
yftm(Rm) = 9#a- This value depends on the nature of the diffusion process and, therefore, 
on the value of the effective diffusion coefficient, which in the considered case can 
be limited to D% < D e ff < D t , where D% and D t are the Bohm and turbulent diffusion 
coefficients, respectively. 

Bohm diffusion The Bohm diffusion coefficient is D# = «b ( c ^B7b^ m /2 e A0> where 
«B is the efficiency parameter, which according to [13] lies in the range 0. 1 — 0.25. 
Putting this value to Eq. (15), combining it with Eq. (13) and solving equation 
9^in(^m) = for R m , one finds 

j£> = 7 x 10' cm x „r T 6 - 2/ " m 5 " 3 (^) . (16, 

where OC0.1 = a/0.1. The plasma temperature at the magneto spheric boundary, T 6 , 
is normalized here according to [27] and expressed in units 10 6 K. Finally, L37 is 
the X-ray luminosity of the pulsar in units 10 37 erg s _1 . 
Turbulent diffusion The turbulent diffusion coefficient can be evaluated as 
D t = a t V s h z , where a t is the efficiency parameter of the turbulence in the slab. 
Putting this value to Eq. (15), combining it with Eq. (14) and solving equation 
9ftin(^m) = tXTa for Rm, one finds 

„(t) 1A 9 1/4 -1/4 1/2 1/4^-1/8 r -1/4 / ^ns \ ^ ,, nx 

J?V ^ 10" cm X < 7] m M 3 m ' T 6 hi (jTj^J • ( 1V ) 

(B) 

The value of Rm is smaller than the Alfven radius of the neutron star evaluated in the 
non-magnetized accretion flow scenario, Ra, by a factor of 5, but remains substantially 
larger than the radius of the neutron star, R ns ~ 10 km for any reasonable value of the 

X-ray luminosity. The condition R$ < R& is satisfied if (a t /r/ m ) hi^ % — 0-02. This 
indicates than the plasma entry into the stellar field in the Magnetically Controlled 
Accretion (MCA) scenario can be explained in terms of diffusion and no additional 
assumptions about instability of the magneto spheric boundary are required. 

5. TORQUE APPLIED TO NEUTRON STARS IN MCA SCENARIO 

For an accreting slowly-rotating neutron star to spin down, the angular velocity of the 
accreting material at the magneto spheric boundary should be smaller that the angular 



velocity of the star itself, i.e. co en < co s [7]. This condition in the case of GX 301-2 is 
satisfied only if 



t, < 0.001 



^orb \ f Rg \ 2 ( Rrn^ 



685s/ l41.5d/ 1 2.5x 10 H cm/ V 10 s 



'cm 



This is two orders of magnitude smaller than the average value of this parameter derived 
in numerical studies of spherical accretion in non-magnetized flow approximation [33]. 
On the other hand, a significant dissipation of angular momentum is expected in the case 
of a magnetized flow (see e.g. [28]). This finding, therefore, favors the MCA scenario in 
this pulsar. 

The absolute value of spin- up torque applied to the neutron star from the accreting 
material, = ^CIo^RqWIz, m m ^ s case > is a factor of 60 smaller than the absolute 
value of the spin-down torque inferred from observations [10], 

K sd = 2nl\v sd \ ~ 6 x 10 32 7 45 ( 10 -i^ s -i ) d y necm > (18) 

and hence, it does not affect the spin behavior of the pulsar during this stage. Here 745 is 
the moment of inertia of the neutron star in units of 10 45 gem 2 . 

The spin-down of an accreting neutron star is associated with interaction between 
the magnetosphere co-rotating with the star and the material at the magneto spheric 
boundary. Since both the Reynolds number and magnetic Reynolds number in the slab 

are very large, this interaction leads to turbulization of the material at the boundary. The 

it] • ? 

spin-down torque associated with this process can be evaluated as K^ d = £^^(0^^, 
where QJt e ff is the mass involved into the turbulent motions in a time unit, and £ t is the 
efficiency parameter. The energy transferred by the turbulent motions in a time unit is, 

correspondingly, E t = co s K^ = e t Wl eS V%(R m ). 

Parameter 9DT e ff can be evaluated by equating the ram pressure of the material in the 
slab with the magnetic pressure due to the dipole field of the neutron star, that is 

u 2 

9*eff = 



7^ /2 (2GM ns )i/2' 

Thus, the spin-down torque applied to the neutron star surrounded by the slab is 

(,) fit H 2 G) s 

Sd ^ 2 (2GM ns )V2- (9) 

(B) 

Putting R m = Rjn into this equation and using parameters of GX 301-2 one finds 

/ O \ 17 /13 / p \ 57/13 
y(0~9 x in33 pry -3/13p-l T 3/13,6/13 -14/13 I B ) ( R ™ ) Hvnerm 

£ sd _2xl0 e t a 01 P 6S5 T 6 L 37 m [b^J [lOk^J d y necm ' 

where is the spin period of the pulsar in the units of 685 s, and 5crsf = 4 x 10 12 G is 
the strength of the surface field of the neutron star measured through observations of the 



Compton Resonant Scattering Feature in its X-ray spectrum ( see [21], and references 
therein). Thus, the spin-down torque applied to the neutron star in GX 301-2 during the 
spin-down trends of the pulsar can be explained in terms of the MCA scenario provided 
the turbulent efficiency is £ t > 0.3, which is in a good agreement with theoretically 
predicted value (see e.g. [23]). 

6. CONCLUSIONS 

The currently used scenario of spherical accretion in HMXBs is solely built around the 
assumption that the accretion flow is non-magnetized. We show this assumption to be 
controversial. Recent results of surface field measurement of O/B-type stars suggest that 
the magnetic field energy density in the material captured by neutron stars in wind-fed 
HMXBs is comparable with (or even higher than) the thermal energy density of the wind 
at the Bondi radius. The accretion flow in this case can be strongly affected by its own 
magnetic field. The magnetic field in the spherical accretion flow is rapidly growing. 
The flow is decelerated as the magnetic pressure due to magnetic field of the flow itself 
reaches its ram pressure. The distance to the deceleration region (Shvartsman radius) 
under certain conditions exceeds the Alfven radius of the neutron star. The character- 
istic time of accretion process inside the Shvartsman radius is limited to the time of 
annihilation of the magnetic field in the accretion flow, which under the conditions of 
interest substantially exceeds the dynamical time as well as the cooling time of the ac- 
cretion flow in this region. The accreting material in this case accumulates around the 
magnetosphere of the neutron star forming a dense non-Keplerian slab. The entry of the 
material from the slab into the magnetosphere is governed by turbulent diffusion. The 
pulsar spin evolution in this Magnetically Controlled Accretion (MCA) scenario is de- 
termined by the process of angular momentum exchange between the star and the slab. 
The rate of angular momentum transfer evaluated in this scenario is in good agreement 
with observations. 
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